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Abstract. - Biased diffusion of two species with conserved dynamics on a 2xL periodic lattice 
is studied via Monte Carlo simulations. In contrast to its simple one-dimensional version on a 
ring, this quasi one- dimensional model surprisingly exhibits phase separation in its steady state, 
which is characterized by one macroscopic particle cluster. We study the order parameter and 
the cluster-size distributions as a function of the system size L, to support the above picture. 



Driven diffusive lattice gases are among the simplest systems exhibiting generic non- 
equilibrium behavior in their steady states 0, a. Even the simplest group of these models 
with purely short-range interactions, produced by the excluded volume constraint, can display 
complex phase diagrams when driven away from equilibrium. The ones with non-trivial phase 
diagram typically involve a breaking of translational invariance (e.g., open boundaries ||), or 
more than one species of particles |j| . The study of multi-species systems is also motivated 
by fast ionic conductors with several mobile species || , water droplets in microemulsions with 
distinct charges M , gel electrophoresis || ^] , and traffic flow |[o| . 

One striking feature of the non-equilibrium steady states of some of these models is that they 
exhibit spontaneous symmetry breaking or phase separation in one dimension (ID), in contrast 
to equilibrium systems with short-range interactions, where no long-range order (LRO) exists 
at finite temperature, i.e. with noise present. Indeed, a simple two species asymmetric exclu- 
sion model with open boundaries display the unusual phenomenon of spontaneous symmetry 
breaking in one dimension |TT| . Also, phase separation, analogous to Bose condensation, 
triggered by defects has been observed in systems with disorder, related to simple traffic flow 
[ [l2| . Recently, a quasi ID model for the sedimentation of colloidal crystals [O, and a simple 
ID "cyclic" three-state lattice gas H], both homogeneous systems with no boundary effects 
(i.e., with ring geometry), have been shown to undergo phase separation. 
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In this Letter we aim to illustrate that the role of dimensionality contributing to long-range 
order and the notion of "lower critical dimension" can be quite subtle in non-equilibrium 
systems. We focus on the difference between the ID and the quasi ID version of a driven 
three-state lattice gas, consisting of holes and two oppositely "charged" species of particles, 
subject to an "electric" field, E. Particle-particle exchanges (PPE) are allowed but occur 
on a much slower time scale than the dominant particle-hole exchanges. In two dimensions 
(2D), the system exhibits a transition from a disordered (homogeneous) to an ordered phase 
characterized by a compact strip of particles spanning the system transverse to the field Q . 
Indeed, even in the absence of PPE, a phase transition exists . In one dimension, the behavior 
of a system without PPE is clearly trivial. With PPE, it was first believed that the ID system 
would mimic those in higher dimensions and display a transition to an ordered state. However, 
an exact solution (with E = oo) was subsequently found |^5| , showing that this system never 
orders (though it exhibits non-trivial cluster size distributions and microscopic shocks 
|fl6|). This behavior seems counter-intuitive, since blockages, which cause the transition to 
an ordered state, should be enhanced in a ID system. On closer examination, blockages are 
found to be "over-enhanced" , in the sense that they occur on microscopic length-scales and 
prevent the macroscopic cluster from growing. As a result, the steady state of the ID system 



is characterized by a typical particle cluster size, which depends only on the microscopies 1 15 
and does not scale with the system size. 

For equilibrium systems with short-range interactions, LRO cannot exist in a ID chain or 
in a (quasi ID) pair of chains. In contrast, the behavior of our model is much more intriguing. 
Using Monte Carlo simulations and simple theoretical methods B, we investigate the steady 
states of a ID system of L sites and the quasi ID case of 2xL. Confirming that there is no 
LRO in the ID case, we find that the 2xL system behaves as an LxL one! In other words, 
a macroscopic cluster forms for sufficiently large E and particle density, with density profiles 
resembling the mean- field ones closely. For maximal effects, we choose E — oo, so that we can 
compare directly with the exactly solvable model in ID fig] . We have also considered a range 
of transverse diffusion rates, as a further attempt to interpolate between the ordering, 2xL 
system and the non-ordering ID case. 

A concise specification of our model consists of a (quasi ID) periodic lattice of 2xL sites. A 
site x = (x, y) can be empty or occupied by either a positive or a negative particle. Associated 
with these occupancies are the standard variables , n^, which assume the value or 1. The 
excluded volume constraint also implies ^x^x = 0' f° r an y x - The external field is chosen so 
that a positive (negative) particle never moves in the —y (+y) direction. At each elementary 
time step, a pair of neighboring sites ("bond") is chosen randomly. For an "x-bond", particle- 
hole pairs are exchanged with rate T±, while PPE occurs with rate 7r_i_- Similarly, for 
"y-bonds", the rates of the allowed exchanges are F\\ and ^\\, respectively. Note that 7 sets 
the time scale for PPE processes. In the simulations presented here, T± = T\\ = 1, 7 = 0.10. 
As indicated, we also studied cases with r±/P|| < 1, since the = limit corresponds to two 
uncoupled ID systems. Our time unit is one Monte Carlo step (MCS), during which 4L bonds 
are chosen. We also restrict ourselves to "neutral" systems (X)x n x = Ex%) a ^ na ^ nnm g> 
i.e., m=j£ X)x( n x + n x) = 0-5- On lattices with L ranging from 10 2 to 10 4 , our runs last 
from 5xl0 5 to 2xl0 6 MCS. Averages, denoted by (• • •), are performed over the time series, 
once the system has settled in a steady state. To ensure that the final state is independent of 
the initial configuration, we started the simulations with both disordered (random) and fully 
ordered configurations. 
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Fig. 1. - (a) Time-series of Q and j + for L = 1000, starting from random initial configuration. 
The inset shows a typical ordered configuration for the same L (bias points to the right); black, 
dark gray, light gray corresponds to +'s, — 's, and holes respectively, (b) Finite-size effect for 
(Q) and ((AQ) 2 ) (inset) in the ordered phase. The solid line in the inset represents a 1/L 
power law. Note the extremely small change in (Q) (less than 0.5% of its saturation value). 
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To characterize the steady state, we define the particle density profile (in y): 
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and an (unnormalized) order parameter, 
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Being the magnitude of the lowest Fourier component of p, Q provides a sensitive measure of 
the density inhomogeneities. For the completely ordered configuration (i.e., p — 0(mL — y); 
being the Heaviside step- function) , Q — sin(m7r)/(Lsin("7r/L)) w 0.318 for an infinite system 
with m = 0.5. On the other hand, (Q) vanishes in the disordered phase, up to finite size effects 
of 0(1/ VL). We measured j + , the current of positive particles. By symmetry, the average 
negative particle current should be just — Finally, we also constructed size distributions 
for both particle and hole clusters by building histograms with respect to their length along 
the field in one fixed "column", say for x = 1. 

For all of our system sizes, we observed that the steady-state configuration is ordered (inset 
of fig. la). Most of the particles "condense" into one macroscopic cluster, while the remainder 
scatters as a small but finite density of "travellers" through the empty region. Due to the 
infinite field, holes cannot enter the macroscopic particle cluster. Note that there will always 
be a finite current of either species (proportional to 7) in the system, so that coarsening does 
take place. The growth of this cluster, when the initial configuration is random, is quite 
interesting. After the first 10-20 MCS, small blockages (particle clusters) form everywhere. 
After this initial phase, a somewhat slower process takes over: particle clusters coarsen in 
time until a single, macroscopic cluster remains in the system. The evolution of Q (fig. la) 
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provides a good picture of this growth process. One important observation is that the value 
of the current is basically unaffected after the initial appearance of small clusters throughout 
the coarsening process (fig. la). This can be understood: the current is mainly controlled by 
PPE within the bulk of the particle clusters. 

Taken alone, the appearance of a single large cluster is not sufficient for us to conclude the 
presence of an ordered state. Indeed, even for the ID case, in which there is a finite typical 
particle cluster size (I* ~ 47 -2 for m = 0.5) we would generally observe a single cluster 
if L <C I*. Although the range of our L's does exceed Aj~ 2 — 400, we probe more deeply by 

analyzing the L dependence of various quantities. The conclusion is that, if jL 1 ^ ~ 1 were the 
criterion for the 2xL system to display ID behavior, then u> ~ 0. In other words, we believe 
that the collective properties of the quasi ID system are distinct from the non-ordering ID 
chain, but fall within the class of the higher dimensional systems. In the remainder of this 
letter, some details of our analyses are provided. 

In fig. lb, we show the average order parameter, (Q), as a function of L. It saturates rapidly; 
even with L = 100, it is within 0.5% of the apparent L = oo limit. Meanwhile, its fluctuations, 
((AQ) 2 ) = (Q 2 ) — (Q) 2 , decay like 1/L (inset fig. lb). The distribution of the sizes of particle 
clusters reveals more of the nature of the steady state. Denoting by P(l, L) the probability 
of finding a cluster of length I in a system of size L, we are not surprised that there are two 
well-separated components, one corresponding to the single macroscopic cluster and the other 
to the small clusters of "travellers". The macroscopic component appears as a peak at Iq. We 
found that Iq w 0.471L, which is a result of m = 0.5 and a clear signal of the large cluster 
scaling linearly with the system size. Meanwhile, 5lo, the standard deviation, grows as VI. 
Focusing on this component alone, we renormalize its integral to unity and rewrite it in terms 
of u = (/ — Iq)/SIq. The result, p(u,L), is well fitted by a Gaussian (fig. 2a). Deviations from 
the Gaussian, seen especially for small L, carry non-trivial information on finite size effects 
which should be investigated further. In the "travellers" component, the distribution decays 
exponentially (oc e~'/ A ), with A < 1 and independent of L. Similarly, the size distribution of 
the hole clusters is also a simple exponential, independent of L apart from finite size effects. 
Both of these exponentials can be understood O] from the observation that the small density 
of "travellers" is roughly homogeneous. Following approaches to the percolation problem Jl8[ , 
we consider another interesting distribution, namely, the probability that a randomly selected 
particle belongs to a cluster of length / (at a given time). Clearly proportional to IP(1,L), 
it will be referred to as the "residence" distribution. To compare distributions from different 
L's, we define w = l/(mL) e (0, 1]. Thus, P(w,L) is the probability (density) that a particle 
"resides" in a cluster of length wmL. Since there are essentially no holes in the particle 
clusters, w is just the fraction of all the particles in such a cluster. In fig. 2b, we see that 
P(w, L) peaks at approximately and at w ~ 0.942. The areas under each component are 
essentially independent of system size, while both sharpen with increasing L. With finite L, 
the former peaks at l/(mL), with area 0.058. Of course, the complement area is under the 
other component. It is not a coincidence that this area (0.942) is identical to the peak position, 
wo, since both are simply related to the fraction of particles in the macroscopic cluster. Finally, 
note that this quantity is also lo/(jnL). 

We also studied the system with varying transverse hopping rates: < Pj_/rii < 1. Our 
results indicate that a crossover to the ID behavior occurs at a very small but finite value of 
this ratio (of the order of 10~ 3 ); only a very careful finite-size analysis could reveal whether 
this occurs when the transverse hopping rate becomes comparable to the inverse relaxation 
time of the ID system. 

Most of our data can be understood by simple mean-field considerations [0 . As we already 
pointed out, the current of one species {e.g., that of the +'s) is controlled by the bulk of 
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Fig. 2. - (a) The distribution of the relative size-fluctuations of the macroscopic particle clusters 
for different system sizes. The solid line is a Gaussian with zero mean and unit standard 
deviation, (b) The scaled residence distributions. System sizes and symbols are the same as 
in (a). The inset gives a closer view of the distribution around wq — 0.942. 



the particle cluster. Although there is a non-trivial charge distribution within the particle 
cluster 0, it is fairly flat (up to finite size effects) around the center of the cluster: thus 
at the simplest level we can approximate the density of both species there to be equal to 
1/2, and « 7/4. Since the steady state is stationary, the current must be constant 

and homogeneous through the system, so that ~ 7/4 outside the particle cluster as 

well. Assuming a homogeneous density of "travellers", m*, outside the macroscopic particle 
cluster, we have f§ 7/4 = (m*/2)[(l - m*) + j(m*/2)], yielding m* = 7/(2-7). For 7 = 0.1, 
m* = 0.0526, which is within 5% of the Monte Carlo value. The mean macroscopic cluster size 
lo can also be deduced, since l + m*(L — l ) = mL. The result, 0.472L, is very close to the 
one obtained from the measured particle cluster size distributions. Also, the order parameter 
for such an ordered configuration is just Q = 0.30, which compares well with the data (fig. 
la). Since the velocity of a tagged particle is approximately 7/2 in the particle cluster and 
(1 — m*) +7(m*/2) outside the cluster, it is easy to show that the typical time it spends 
diffusing through the macroscopic cluster is a fraction l / (mL) of the total time needed to 
travel through the system once. Thus, for our parameters, a randomly selected particle is 
found with probability 0.944 in the large particle cluster, which reproduces the location of the 
second peak of the distribution P(w, L) quite well. 

In summary, we have studied the non-equilibrium steady state of a driven, quasi one- 
dimensional, three-state stochastic lattice gas, using simulations and simple mean-field argu- 
ments. We found that, for the parameter regime studied, the properties of a coupled pair 
of chains of length L are drastically different from those of the strictly ID system, i.e., the 
single chain. With as little as 1% of cross-chain (transverse) moves, the 2xL model already 
develops a single macroscopic blockage, similar to the 2D (LxL) system. Moreover, the 
finite-size analysis of the ordered phase, up to L = 10 4 , shows no indication of a crossover to 
ID behavior. Thus, we conjecture that LRO survives the L — > 00 limit in our simple model 
(fixed 7 and T± = If this conjecture proves to be true, our model would pose a stark 

contrast to, e.g., the MxL equilibrium Ising model where LRO does not exist in the fixed 
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M, L — ► oo limit. Instead, the strictly one-dimensional system is a "singular limit", in much 
the same way that entirely different behaviors are found, obviously, in IxL vs. 2xL systems 
of biased diffusion of two species without particle-particle exchange in finite field. The key 
difference between the strictly and the quasi ID systems appears to be the ease with which 
particles can pass one another, leading to the formation of a macroscopic cluster. Work is in 
progress to reveal further differences between the ID and the quasi ID systems, e.g., details 
of the coarsening process leading to LRO (TtJ . 

Note Added. - After submitting our manuscript we became aware of the work by Arndt 
et al. [ fl9| |. They studied a ID model on a ring, similar to ours, but with different exchange 
rates. 

#*# 

Stimulating discussions with S. Sandow, C. Godreche and D. Mukamel are gratefully ac- 
knowledged. Special thanks to Z. Toroczkai for his comments on the manuscript. We thank 
E. Ben-Nairn and G. Lythe for their support through the DOE at the CNLS Nonequilibrium 
Dynamics Workshop (April 1998) where this work was first presented. G. Korniss also wishes 
to thank P. A. Rikvold and M. A. Novotny for their support during the preparation of this work 
and their critical reading of the manuscript. This research at Virginia Tech was supported in 
part by NSF Grant No. DMR-9727574. Work at the Florida State University was supported 
by SCRI under the US DOE Contract No. DE-FC05-85ER25000, and by NSF Grant No. 
DMR-9520325. 



REFERENCES 

Katz S., Lebowitz J. L. and Spohn H., Phys. Rev. B., 28 (1983) 1655; J. Stat. Phys., 34 
(1984) 497. 

Schmittmann B. ZlA R. K. P., in Phase Transitions and Critical Phenomena, edited by Domb 
C. and Lebowitz J. L., Vol. 17 (Academic Press, N.Y.) 1995. 

Krug J., Phys. Rev. Lett, 67 (1991) 1882; Schutz G. and Domany E., J. Stat. Phys., 72 (1993) 
277; Derrida B., Evans M. R., Hakim V. and Pasquier V., J. Phys. A, 26 (1993) 1493. 
Schmittmann B., Hwang K. and Zia R.K.P., Europhys. Lett., 19 (1992) 19. 
Korniss C, Schmittmann B. and Zia R. K. P., Europhys. Lett., 32 (1995) 49; J. Stat. Phys., 
86 (1997) 721. 

See e.g., Chandra S., Superionic Solids. Principles and Applications, (North Holland, Amster- 
dam) 1981. 

Aertsens M. and Naudts J., J. Stat. Phys., 62 (1990) 609. 

Rubinstein M., Phys. Rev. Lett, 59 (1987) 1946; Duke T. A. J., Phys. Rev. Lett, 62 (1989) 
2877; Schnidman Y., in Mathematics in Industrial Problems IV, edited by Friedman A., 
(Springer, Berlin) 1991; Widom B., Viovy J. L. and Desfontaines A. D., J. Phys I, 1 (1991) 
1759. 

Alon U. and Mukamel D., Phys. Rev. E, 55 (1997) 1783. 

Biham O., Middleton A. A. and Levine D., Phys. Rev. A, 46 (1992) R6128; Leung K.-t., 
Phys. Rev. Lett, 73 (1994) 2386. 

Evans M. R., Foster D. P., Godreche C. and Mukamel D., Phys. Rev. Lett, 78 (1995) 
208; J. Stat. Phys., 80 (1995) 69. 

Evans M. R., Europhys. Lett, 36 (1996) 13; Krug j. and Ferrari P. A., J. Phys. A., 29 
(1996) L213. 

Lahiri R. and RAMASWAMY S., Phys. Rev. Lett., 79 (1997) 1150. 

Evans M. R., Kaf ri Y., Koduvely H. M. and Mukamel D., Phys. Rev. Lett, 80 (1998) 425; 



cond-mat/980225E Preprint, 1998; see also Arndt P. F., Heinzel T. and Rittenberg V., J. 
Phys. A., Al (1998) L45. 



G. KORNISS et al: ■ LONG-RANGE ORDER IN A QUASI ONE-DIMENSIONAL ETC. 



7 



[15] C. Godreche and S. Sandow, (private communications, to be submitted). 

[16] Derrida B., Janowsky S. A., Lebowitz J. L. and Speer E. R., Europhys. Lett., 22 (1993) 

651; J. Stat. Phys., 73 (1993) 813. 
[17] Mazilu L, Korniss G. and Schmittmann B., (in preparation). 

[18] See e.g., Stauffer D., Introduction to Percolation Theory, (Taylor & Francis, London and 

Philadelphia) 1985. 
[19] Arndt P. F., Heinzel T. and Rittenberg V. 



:ond-mat/9809123| Preprint, 1998. 



